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Abstract 


We have completed the 0{as) QCD corrections to exclusive heavy quark-antiquark distri¬ 
butions in deep inelastic electroproduction and present here the differential distributions in 
the masses of charm-anticharm and bottom-antibottom pairs at the electron-proton collider 
HERA. 
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Order as QCD corrections to the structure functions for single particle inclusive deep 
inelastic electro-production of heavy quarks were recently published in |l|. The reaction 
e~(h) + P{p) —^ + Q{pi) + Q{P 2 ) -I- X is dominated by the virtual photon mediated 

reaction when —= —(/i — I 2 Y ^ ^'zi the heavy quark differential production cross 
section can be written as 

+ (1 - yf}F 2 {.x, rri^) - y^Fiix, Q^, m^)] , (1) 

after integration over the azimuthal angle between the plane containing the incoming and 
outgoing electron and the plane containing the incoming proton and outgoing heavy quark. 
The square of the center of momentum energy of the electron-proton system is denoted 
by S, and the variables x and y are dehned as x = /2p ■ q and y = p ■ q/p ■ h with 

—= xyS. The heavy quark structure functions F 2 {x,Q‘^,m‘^) and Fl^x^Q"^ 
are functions of the heavy quark mass m. We assume that the heavy quark production 
is extrinsic so that F 2 and Fl can be calculated from an analysis of the virtual photon 
induced reaction 7*(g) -|- P{p) —>■ Q{pi) + Q{p 2 ) + X and its corresponding parton analogue 
7 *(g) -|- ai(/ci) —>• Q{pi) -l- Q{p 2 ) + ^ 2 (^ 2 ), where Oi and 02 are zero-mass gluons g or light 
mass (anti) quarks (g) g as opposed to the massive (anti) quarks {Q) Q. The result is that 
the structure functions can be obtained from the partonic results via the formula 


Fk{x,Q ,m ) = 


Q^Osip^) P 2 




^min 
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where {k = 2, L). The lower boundary on the integration is given by ^min = x{4np + Q'^)/. 
Further /j(^, ,{i = g, q, q) denote the parton momentum distributions in the proton, and p 

stands for the mass factorization scale which has been put equal to the renormalization scale 
in the running coupling constant. Finally, and , (z = g ,q ,q',l = 0,1), and d^j}] and 

, {i = q ,q) are scale independent parton coefficient functions which were hrst calculated 
in [|I|. In Eq. (2) we made a distinction between the coefficient functions with respect to 
their origin. The coefficient functions indicated by and c^\- originate from the partonic 
subprocesses where the virtual photon is coupled to the heavy quark hence the factor of ejj. 
The quantity d^j}] comes from the subprocess where the virtual photon interacts with the 

light quark so it is proportional to e^. The quantity comes from the interference between 
the above processes and hence has a factor enGi with all charges in units of e. Note that 
terms proportional to cnGi appear in the photon-parton differential distributions even though 
they integrate to zero in the total partonic cross section. Furthermore we have isolated the 
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factorization scale dependent term containing ln(/i^/m^). The functions multiplied by this 
term, which are indicated by a bar, are called mass factorization parts. Note than Eq. (2) 
only holds for > 0. In the photo-production limit there are additional terms involving 
the parton densities in the photon. 

The previous treatment of the 0{as) corrections yielded results for the inclusive distribu¬ 
tions for heavy quarks in the virtual photon induced reaction 7 * -|-P(p) —>• Q{pi){Q{p 2 )) + ^, 
i.e., the differentials dFk{x,Q‘^,m‘^,pt)/dpt and dF^^x, , rri^, y)/dy. The corrections to the 
heavy quark inclusive pt and y distributions at hxed points in and x were published in 
P]. Event rates for regions of the x and plane have been presented in [^. 

In this paper we report on the results of a calculation of the 0{as) corrections for heavy 
quark exclusive distributions at fixed and x, which allows us to study all correlations 
between the outgoing particles in the virtual photon initiated reaction 7 *(g) -|- P{p) —>■ 
Q{Pi)fQ{P 2 )+X (^ 2 ) with X = 0 or 1 jet (massless parton). This information is of immediate 
interest to the experimenters working with the HI and ZEUS collaborations at the electron- 
proton collider HERA. We therefore present the effects of the QCD corrections to invariant 
mass distributions of a heavy quark-antiquark pair, for 8.5(GeV/c)^ < < 50(GeV/c)^ 

and 4.2 x 10“^ < x < 2.7 x 10“^. The (L, 2 ) photon components are treated separately. As 
input we use the latest GTEQ parton densities 0, which £t the newly released HERA data 

i- 


Our new analysis of exclusive heavy quark deep inelastic electro-production at HERA 
extends the existing studies of inclusive QGD corrections in the virtual photon channel [Q, 
0], 0, inclusive QGD corrections in the real photon = 0) channel 0, [0, and exclusive 
QGD corrections in the real photon channel |^, allowing for an extensive comparison with 
present and future experimental data. Heavy quark electro-production is expected to play an 
important role in the determination of the gluon distribution function in the proton at small 
X. A knowledge of the production cross sections and distributions for charm and bottom 
quarks is also relevant in the study of the GKM matrix elements through the rare decays of 
D— and B— mesons and the analysis of DD and BB mixing |p. 

In our exclusive computation we use the same techniques as the authors of [§ for com¬ 
puting heavy-quark correlations in photo-production and hadro-production. These are based 
on the replacement of divergent terms in the squared matrix elements by generalized plus 
distributions. The divergent terms appear when the propagators diverge in regions of phase 
space where the outgoing parton is soft and/or collinear to the propagating particle. The 
replacement of the divergent terms by generalized plus distributions allows one to isolate 
the soft and collinear poles within the framework of dimensional regularization, without 
having to calculate all the phase space integrals in a spacetime dimension n 7 ^ 4 as usually 
required in a traditional inclusive computation. The resulting expressions for the squared 
matrix elements appear in a factorized form where poles in n — 4 multiply splitting functions 
and lower order squared matrix elements. The cancelation of singularities is then performed 
using the factorization theorem |M. Since the hnal result is in four-dimensional space time. 


we can compute all relevant phase space integrations using standard Monte Garlo integra¬ 
tion techniques and produce histograms for exclusive, semi-inclusive, or inclusive quantities 
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related to any of the outgoing particles. We therefore have a new calculation of the scale 
independent coefficient functions and We checked the rj = s/Am? — 1 

and ^ dependence of the scale independent coefficient functions against the results 

in |I|. The analogous results for the inclusive distributions, dFk/dpt and dFk/dy, were also 
checked against the results in |Q. Additional distributions and correlations along with details 


of the calculation will be presented in a more complete article |11 


Folding the parton densities in the proton with our new scale independent coefficient 
functions as dictated by Eq. (2), we present results for the differential structure functions in 
the invariant mass of the heavy quark-antiquark pair, which we will call M. Thus we give 
plots of dF 2 {x,Q‘^,m‘^,M)/dM and dFL{x,Q‘^,m‘^, M)/dM at fixed x and Q^. We use m = 
rric = 1.5 GeV/c^ for charm production and m = mi, = 4.75 GeV/c^ for bottom production. 
We choose the factorization (renormalization) scale as + 4(m^ -|- {Pt^ + Pt-)^/4) for 

charm production and p? = + Pt^)‘^/A for bottom production. As mentioned 

earlier we use the GTEQ3M parton densities Q in the MS scheme and the two loop a* with 
A 4 = 0.239 GeV for charm and A 5 = 0.158 GeV for bottom. 

Tables 1 and 2 show the variation of F 2 {x,Q‘^,m‘^) and Fi{x,Q'^,m?) with the renor¬ 
malization scale for charm production, and Tables 3 and 4 show the variation for bottom 
production. For charm production typical variations from the central value are less than 
15% and for bottom production they are less than 6 %. 

Figures 1 and 2 show the distributions dF 2 {x, Q^, m^, M)/dM and dFi^x, m^, M)/dM 
for charm production sX x = 8.5 x 10“^ while varying . Figure 3 and 4 shows the 
distributions dF 2 {x,Q‘^,ml, M)/dM and dFL{x,Q‘^,ml, M)/dM for various values of x at 

= 12 (GeV/c)^. For charm production, the Born result multiplied by a constant factor 
gives good agreement with the complete 0{a1) result for the 2 projection, while the L pro¬ 
jection is not well reproduced for a constant multiplicative factor for large invariant masses. 
Figures 5 and 6 show the distributions dF 2 {x,Q‘^,ml, M)/dM and dFL{x,Q‘^,ml, M)/dM 
for bottom production at x = 8.5 x 10“^ while varying Q^. Figure 7 and 8 shows the distribu¬ 
tions dF 2 {x,Q‘^,ml,M)/dM and dFL{x,Q'^,ml, M)/dM for bottom production for various 
values of x at = 12 (GeV/c)^. For bottom production, the Born result multiplied by a 
constant factor gives quite good agreement with the complete 0(a‘^) result for both the 2 
and L projections. 

To conclude, we repeat that the invariant mass distributions are reasonably well repre¬ 
sented by taking the Born result times a multiplicative factor. The agreement is excellent 
for bottom production and not quite so good for charm production. 
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Table 1 


X 


/i = /io/2 

fi = fio 

fi 2/ig 

8.5 xlO-4 

8.5 




8.5 xlQ-^ 

12 

0.97 xlO-i 

1.10 xlO-1 

1.20 xlO-^ 

8.5 xlO-4 

25 

1.57 xlO-^ 

1.79 xlO-i 

1.94 xlO-i 

8.5 xlQ-^ 

50 

2.25 xlO-^ 

2.52 xlO-i 

2.71 xlO-i 






4.2 xlO"^ 

12 


1.38 xlO”^ 


8.5 xlO-4 

12 

0.97 xlO-i 

1.10 xlO-i 

1.20 xlO-i 


12 





12 





Table 1. Variation of F 2 with = Q‘^ + 4(m^ + {Pt^ + Pt-)^/4) for various x and Q‘^ values. 
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Table 2 


X 


/i = /io/2 

fi = fio 

fi 2/ig 

8.5 xlO-4 

8.5 




8.5 xlQ-^ 

12 




8.5 xlO-4 

25 




8.5 xlQ-^ 

50 

QjjoQljBjl 








4.2 xlO"^ 

12 




8.5 xlO-4 

12 





12 

ISJjQJjBH 




12 


li^gQjjyil 



Table 2. Variation of with /ig = + 4(m^ + + Pi-)^/4) for various x and Q'^ values. 
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Table 3 


X 


/i = /io/2 

fi = fio 

fi 2/ig 

8.5 xlO-4 

8.5 



IgjOQljBjll 

8.5 xlQ-^ 

12 


2.45 xlO-3 


8.5 xlO-4 

25 


6.34 xlO-^ 


8.5 xlQ-^ 

50 









4.2 xlO"^ 

12 




8.5 xlO-4 

12 





12 





12 





Table 3. Variation of F 2 with = Q'^ + ml + {Ptf^ + /A for various x and values. 
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Table 4 


X 


/i = /io/2 

fi = fio 

fl ‘^flQ 

8.5 xlO-4 

8.5 




8.5 xlQ-^ 

12 

1.09 xlO-^ 


1.00 xio-^ 

8.5 xlO-4 

25 

4.65 xlO-^ 


4.50 xlO-^ 

8.5 xlQ-^ 

50 

1.57 xl0-‘^ 








4.2 xlO"^ 

12 

1.44 xlO"^ 


1.39 xlO"^ 

8.5 xlO-4 

12 

1.09 xlO-^ 

|2jgQ|jg§ 

1.00 xlO-4 


12 


lO^^QJjgg 



12 


^^joQjjgl 



Table 4. Variation of Fl with /i^ = + {Pt^ + Pt^Y/A for various x and values. 





































Figure Captions 


Fig.l. The distributions dF 2 {x, Q^, M)/dM for charm production at fixed x = 8.5 x 10 ^ 
and = 8.5 (solid line), 12 (dotted line), 25 (short dashed line), 50 (long dashed line) 
all in units of (GeV/c)^. Histograms are complete 0(ag) result. For comparison we 
show the Born result times a multiplicative factor of 1.3 (empty box), 1.3 (solid box), 

1.2 (empty circle), 1.2 (solid circle) for the various values, respectively. 

Fig.2. The distributions (iF/;,(a;, TUc) charm production at hxed a; = 8.5x10 ^ 

while varying Q^. Notation is that of Figure 1 for the 0(ag) result. For comparison 
we show the Born result times a multiplicative factor of 1.9 (empty box), 1.9 (solid 
box), 1.6 (empty circle), 1.6 (solid circle) for the various values, respectively. 

Fig.3 The distributions dF 2 {x, Q^, ml, M)/dM for charm production at hxed = 12 (GeV/c)^ 
with X = 4.2 X 10“"^ (solid line), 8.5 x 10“"^ (dotted line), 1.6 x 10“^ (short dashed line), 

2.7 X 10“^ (long dashed line). For comparison we show the Born result times a multi¬ 
plicative factor of 1.3 (empty box), 1.3 (solid box), 1.3 (empty circle), 1.4 (solid circle) 
for the various x values, respectively. 

Fig.4 The distributions dFii^x, Q^, ml, M)/dM for charm production at fixed = 12 (GeV/c)^ 
while varying x. Notation is that of Figure 3 for the 0(al) result. For comparison we 
show the Born result times a multiplicative factor of 1.9 (empty box), 1.9 (solid box), 

1.9 (empty circle), 1.5 (solid circle) for the various x values, respectively. 

Fig.5. The distributions dF 2 {x,Q‘^,ml,M)/dM for bottom production at hxed x = 8.5 x 
10“"^ while varying Q^. Notation is that of Figure 1 for the 0(al) result but all points 
are 1.3 x Born result for the various values, respectively. 

Fig.6 . The distributions dFL{x,Q‘^,ml,M)/dM for bottom production at hxed a; = 8.5 x 
10“^ while varying Q^. Notation is that of Figure 1 for the 0(al) result but all points 
are 1.3 x Born result for the various values, respectively. 

Fig.7 The distributions dF 2 {x,Q‘^,ml, M)/dM for bottom production while varying x at 
hxed = 12 (GeV/c)^. Notation is that of Figure 2 for the 0{al) result but all points 
are 1.3 x Born result for the various x values, respectively. 

Fig.8 The distributions dFL{x,Q‘^,ml, M)/dM for bottom production while varying x at 
hxed = 12 (GeV/c)^. Notation is that of Figure 2 for the 0{al) result but all points 
are 1.3 x Born result for the various x values, respectively. 


9 



References 

[1] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B392 (1993) 
162. 

[2] E. Laenen, S. Riemersma, J. Smith and W. L. van Neerven, Nucl. Phys. B392 (1993) 
229. 

[3] S. Riemersma, J. Smith and W. L. van Neerven, Rates for Inclusive Deep-Inelastic 
Electroproduction of Charm Quarks at HERA, SMU-HEP 94-25, to appear in Physics 
Letters B. 

[4] H. L. Lai, J. Botts, J. Huston, J. G. Morfin, J. P. Owens, J. W. Qiu, W. K. Tung and 
H. Weerts, MSU-HEP-41024 (Oct. 1994). 

[5] HI Collaboration, 1. Abt, et ah, Nucl. Phys. B407 (1993) 515; ZEUS Collaboration, M. 
Derrick et ah, Phys. Lett. B316 (1993) 412. 

[6] R.K. Ellis and P. Nason, Nucl. Phys. B312 (1989) 551; S. Frixione, M. L. Mangano, P. 
Nason and G. Ridolh, CERN-TH 7527/94. 

[7] J. Smith and W. L. van Neerven, Nucl. Phys. B374 (1992) 36; W. L. van Neerven, 
Nucl. Phys. B29A (Proc. Suppl.) (1992) 199. 

[8] S. Frixione, M. L. Mangano, P. Nason, and G. Ridolh, Nucl. Phys. B412 (1994) 225; 
M. L. Mangano, P. Nason, and G. Ridolh, Nucl. Phys. B373 (1992) 295. 

[9] G. A. Schuler, Nucl. Phys. B299 (1988) 21; A. Ah, G. Ingelman, G.A. Schuler, F. 
Barreiro, M.A. Garcia, J.F. de Troconiz, R.A. Eichler and Z. Kunszt, Heavy quark 
physics at HERA, Proc. DESY workshop on HERA physics, DESY, Hamburg(1988), 
DESY 88-119; A. Ali and D. Wyler in Physics at HERA, vol 2 p. 669, Proceedings of 
the Workshop Hamberg, Oct 29 (1991) Eds. W. Buchmiiller and G. Ingelman; A. Ali, 
DESY 93-105. 

[10] J. C. Collins, D. E. Soper, and G. Sterman in Perturbative QCD ed. A. H. Mueller ( 
World Scientihc, 1989). 

[11] B. W. Harris and J. Smith, in preparation. 


10 



dFj/dM 














•lO ■ 























